Simulating many-body lattice systems on a single nano-mechanical resonator 
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We show that lattice systems, such as the Bose-Hubbard model, can be simulated on a single 
nano- or micro-mechanical resonator, by exploiting its many modes. The on-site Hamiltonians 
are engineered by coupling the mechanical modes to the modes of a pair of optical or striplinc 
resonators, and the connections between the lattice sites are engineered in a similar way. The lattice 
network structure is encoded in the frequency components of the fields driving the resonators. This 
three-resonator configuration also allows universal quantum computing on the nano-resonator. 
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Current methods to simulate many-body lattice sys- 
tems using engineered mcsoscopic systems achieve this 
by fabricating a lattice of the latter |l|~[5[ ■ Here we show 
how to create an effective many-body lattice using a to- 
tal of three mesoscopic resonators. The bodies, or lattice 
sites, are the multiple modes of a nano-mechanical res- 
onator 0, 0| ■ The nonlinear Hamiltonian for each mode 
is engineered by coupling the nano-resonator to a pair of 
optical or superconducting resonators: pairs of the aux- 
iliary modes create the nonlincarity for each mechanical 
mode. The interactions (links) between lattice sites are 
also obtained by a coupling to the modes of an auxiliary 
resonator, which may be one of the previous pair, or a 
third auxiliary. Each lattice link is created by one mode 
of this auxiliary, the structure of the lattice being deter- 
mined by the field that drive it: each frequency compo- 
nent in the driving field selects which nano-mechanical 
mode (or "nano- mode") is coupled to which auxiliary 
mode, and thus which nano-modes arc coupled together. 
This configuration dramatically reduces the number of 
mesoscopic systems required to simulate a many-body 
system, and also means that the lattice network can be 
"reprogrammcd" by changing the driving field. 

Here we focus on constructing the Bose-Hubbard (BH) 
model [H, but we also discuss other lattice Hamiltonians 
that could be engineered in a similar way, including spin 
lattices. The BH model is especially interesting because 
of it's simplicity, and the fact that it possesses both Mott- 
insulator and superfluid phases. It consists of a chain of 
oscillators, with nearest-neighbor linear coupling, each 
with a Kerr nonlincarity. The Hamiltonian for the BH 
model is 

#BH = £( a n a «) 2 + C(OnO„+i + <4i a n+l), (1) 
n n 

where {a n } arc the mode annihilation operators. 

In recent years it has been shown that the modes of 
a doubly-clamped nano-resonator are intrinsically cou- 
pled via a nonlinear interaction, although very weakly Q . 
While this coupling is intriguing, as it can be enhanced 
by driving and exploited in various ways [ToL [111 ] , it is 



still too weak in present systems for our purposes, since 
we wish to work down to the singlc-phonon level. To 
engineer a desired lattice Hamiltonian connecting the 
modes of a nano-resonator, we must therefore couple the 
resonator to another system. Nonlinearities for a sin- 
gle nano-mode, and coupling between nano-modes can 
be engineered by coupling to a non-linear system (see, 
e.g. 12- 16l|) such as a transmon qubit 17, HI- The prob- 
lem with this technique is that if we couple all the modes 
to a single qubit we cannot avoid coupling all the modes 
together simultaneously. While creating an all-to-all net- 
work is an interesting prospect in itself, the above pro- 
cedure can only couple the nano-modes in a lattice con- 
figuration by using many individually fabricated qubits. 
We instead create a lattice by using the multiple (linear) 
modes of a single auxiliary resonator as the set of auxil- 
iary systems to create the lattice links. Then a method 
discovered recently by Ludwig et al. [3] and Stannigel et 
al. [2fj|, building on recent experimental and theoretical 
work [21|, allows us to use these linear modes to cre- 
ate the required nonlincarity for each of the nano-modes 
(with the addition of a second auxiliary resonator) . The 
method works by coupling a single nano-mode to two 
auxiliary modes with the same frequency, while linearly 
mixing the auxiliary modes. There are many ways to 
linearly mix modes, and examples include optical modes 
coupled by a polarization rotator; optical toroidal res- 
onators coupled via an evanescent field [2 if: or super- 
conducting resonators coupled electrically [22j. As we 
elucidate, each pair of mixed modes acts effectively like 
a single qubit, and it is the coupling to these qubits that 
allows us to create the nonlinearities (l5| . 

For notational simplicity, all our Hamiltonians will be 
scaled by 1/h, so that they have units of inverse seconds 
(or equivalently rad/s). With this scaling, all the rate- 
constants that appear in the Hamiltonian are also the 
rate-constants that appear in the resulting equations of 
motion. We therefore write the Hamiltonian for a nano- 



,t, 



where 



mechanical resonator as H m = ^ r 
and a n are the respective frequency and annihilation op- 
erator for each mode. The frequency spectrum is usually 
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highly complex, and some modes will have much higher 
quality factors than others. This is not crucial for our 
purposes, however, as we will be able to select the modes 
we use by choosing the driving frequencies. 

We must also take into account that all our resonators 
are subject to damping. We describe this damping using 
the standard Markovian master equation in the Lind- 
blad form [23|]. Each source of damping or decoherence 
is described by including a term kT>(A)p in the equa- 
tion of motion for the density matrix p, where T>(A)p = 
ApA^ - (l/2)(A*Ap + pA^A) 0. The rate constant k 
is called the damping rate, and the dimensionless oper- 
ator A is called the Lindblad operator. The Lindblad 
operator for the damping of each mode is simply the an- 
nihilation operator for that mode, and their respective 
damping rates will be defined as needed. 

We first consider how to couple selected mechanical 
modes together using the modes of an auxiliary (optical 
or microwave) resonator. If we denote the frequencies and 
annihilation operators of the auxiliary modes by Qj and 
bj, then the auxiliary Hamiltonian is H aux = y"V flMbi, 
and the basic coupling to the nano-resonator is [2 51 ] 



(2) 
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in which every mechanical mode is coupled to every aux- 
iliary mode by some rate-constant g n j . We now create an 
effective linear interaction by using the standard method 
employed in sideband cooling driving the auxiliary 
with a coherent field — but we modulate this drive so 
that it contains N frequencies v^, k = 1, . . . , N. The aux- 
iliary field will settle down to some steady-state, given by 
a coherent part plus quantum fluctuations. Denoting the 
coherent part by (3(t) = J2k Pk exp(— iukt), we remove it 
by transforming to the "displacement picture" [24| . The 
interaction Hamiltonian becomes 



H[ = 9m [b] + P* (*)] [bj + /?(*)] (a„ + 4). (3) 



We now move to the interaction picture with respect to 
the free Hamiltonians of the two resonators, and per- 
form the rotating- wave approximation (RWA) (26j . It is 
this that selects which modes will interact, by eliminating 
pairs that do not satisfy the three-way resonance condi- 
tion uj n + flj — Ufc = between the modes and the driving 
frequencies. The result is H lnt i =J2k ^k, with 



Hk = ^Y^9k\Pk\{akb\ + b k a\), 



(4) 



We now select the driving frequencies so that two me- 
chanical modes, a^\a^ are coupled to the same auxil- 
iary mode bk, and then adiabatically eliminate the aux- 
iliary. To perform this adiabatic elimination we first de- 
tune the nano-modes from resonance with the auxiliary 
by an amount Acj. The adiabatic elimination requires 
that this detuning and the damping rate of the auxiliary, 
k, are much faster than the dynamics of the mechanical 
modes, which in our interaction picture are given purely 
by the coupling rate gk- We can now set bk = 0, and this 
"slaves" the state of the auxiliary to that of the nano- 
modes so that bk = gkPki^ + o^')/{k — iAuj) + h(t). 
Here h(t) is quantum input noise that comes from the 
damping j27|. After applying the RWA once more, we 
obtain the simple linear coupling 
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Since the auxiliary mode is damped, when eliminated 
it also induces an effective damping in the mechan- 
ical modes. The effective damping rate is = 
{gkPk) 2 K(Auj 2 + k 2 ), and the Lindblad operator is Ak = 



j_ a ( 2 > 



Note that since we can couple any number of 
nano-modes to a single auxiliary mode, we can use the 
above technique to couple any subset of lattice sites in 
an all-to-all network if we wish. 

Implementing the above coupling procedure gives us 
a many-body lattice of harmonic oscillators (or bosons). 
To make the lattice really interesting we need to add a 
nonlinearity to each oscillator. As mentioned above, we 
do this using the technique in [l9|, [2(| . We couple each 
mechanical mode to a pair of linearly coupled auxiliary 
modes of the same frequency, whose annihilation opera- 
tors we denote by c& and dk- The coupling Hamiltonian 
is the same as that in Eq.([2]), except that the auxiliary 
modes now come in pairs, and the Hamiltonian contains 
linear mixing of these pairs. The interactions between all 
the modes are therefore given by 

H 2 =5^/nj(^c J -+dt^)(an + aJ l ) + aj(2,-4 + ^)' ( 6 ) 

nj 

where Sj is the rate of mixing for each linear pair. To 
begin we rewrite the full Hamiltonian in terms of the 
transformed modes cj = (dj + dj)/y/2 and dj = (dj — 



dj)j\/2 (these are referred to as "supermodes" [21(), and 
the result is 



H = bj n a) n a n + ^2 [(% + s j) c } c j + (% _ s j)d]d. 
n j 



(7) 



where each k labels a pair of modes (ak,bk). Note that 



the two additional terms multiplying [a n - 



Eq.Q, 



namely b\bk + \/3k\ 2 , have averaged out under the RWA 
applied with respect to the mechanical frequency. 



This transformation to the new modes c and d is reveal- 
ing. It shows us that we have two auxiliary modes now 
at two different frequencies (flj ± Sj ) , and that the non- 
linear coupling to the nano-modes flips photon between 
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these two modes at a rate proportional to (a n + a^) pjj ]. 
Further, the interaction has a clear resonance condition: 
in the interaction picture the term f n jdj(y-a n and its 
Hermitian conjugate will be constant if 2s j = uj n . By 
varying the mode mixing rate Sj across the modes, a 
single auxiliary mode can be selected to interact with 
a single mechanical mode. If one is able to modulate 
the non- linear mechanical/electrical coupling rate f n j at 
a frequency 8 n j, then the resonance condition becomes 
2s j = uj n + 8 n j, providing a second means of tuning the 
mode selection. 

Now focus on the interaction between a single nano- 
mode and auxiliary mode-pair, where we have detuned 
the nano-mode from the resonance by A so that 2sj = 
io n + A. Assuming that all other nano-modes are much 
further from resonance with the mode-pair than the one 
we have selected, applying the RWA to the Hamiltonian 
above, and moving into the interaction picture with re- 
spect to the nano-mode frequency u; n , we have 

H' = f2(c f c + Sd) + A(c t c - d f d) + f[adJ + {ad)% 

We note that since e'e + d^d commutes with all the other 
operators in H' , the high frequency f2 has no effect on 
the interaction dynamics, and thus on the perturbation 
analysis that we will shortly use. But it will be impor- 
tant when we adiabatically eliminate the cavity mode 
c. The Hamiltonian H' is especially interesting because 



the operators dc 1 and 



d^d have the commutation 



relations of spin operators — they are the Schwinger 
representation of angular momentum, with identification 
a z = c)c — d^d and <r+ = dc^ . The linear coupling of the 
modes has converted the resonator's basic nonlinear in- 
teraction with the linear modes into a linear interaction 
with a nonlinear system. Because the nano-resonator 
mode is now effectively coupled to a spin- 1/2 system, via 
a spin operator that does not commute with the spin 
Hamiltonian (in this case cr z ), it is possible at least in 
theory to en gine er a range of non-linear Hamiltonians 
for the mode [la ]. 

The Kerr (or x ) nonlinearity, can now be created 
in at least two ways. We can make A sufficiently large 
(A ^> cj m , /) that the interaction is a perturbation for the 
qubit. If we use time- independent perturbation theory 
to diagonalize the Hamiltonian in powers of e = (//A), 
the even powers give the effective interaction a m (a^a) 2m , 
m = 1,2,..., where a m is some Pauli spin-operator for 
the qubit. The forth-order term (m = 2) is the Kerr 
term. However, to create the Kerr in this way requires 
that we prepare the qubit in the correct eigenstate of 
a m , and this is not so easy for our effective optical qubit. 
The other method, used in [ljj [2(|, is to expand the 
perturbative series to second order, giving the interac- 
tion H e g = <Ti(aJa), and then adiabatically eliminate the 
qubit. This is essentially the same adiabatic elimination 
procedure we used above. In this case the qubit opera- 
tor o\ becomes effectively (a^a) , so that H e s becomes the 



Kerr Hamiltonian. 

To perform the perturbative diagonalization to second- 
order, we transform the Hamiltonian using U = 
cxp[—i(adc ! — aWc)/2], The resulting effective Hamil- 
tonian commutes with both e'e and d^d, so it has the 
same eigenstates as the latter. Leaving mode d in the 
vacuum state, all appearances of d^d disappear, and we 
are left with the very simple non-linear Hamiltonian 



H" =ficV 



(8) 



Due to the damping of the auxiliary modes at the rate k, 
the transformation U also creates an effective damping 
for the resonator mode at rate 7' = e 2 k/4, with Lindblad 
operator Ly = ca. 

We now drive mode c with a field a, shift to the 
displacement picture to create the effective interaction 
a(c + c ! )(a'a) (the same procedure that gives Eq.([3|) 
above), and adiabatically eliminate the auxiliary mode 
c. It is now useful to remember the following general 
rule: given a mode 6 that interacts with another system 
via (b + b')K, where b has a fast oscillation rate f2 and 
a fast damping rate K, then to second-order in the adia- 
batic elimination, f2 produces the effective Hamiltonian 
K 2 (in a manner similar to perturbation theory), and 
the damping k induces an effective damping (or continu- 
ous measurement) via the Lindblad operator K . For the 
above system the result is the Hamiltonian [2(J 



H 



Kerr 



(efa/A) 2 n 

n 2 + n 2 



(a f a) 2 , 



(9) 



with decoherence at rate T = (e f 'a j '4) 2 n / '(ft 2 + n 2 ) via 
the Lindblad operator L = a T a. 

The two effective Hamiltonians we have described, 
Eqs.([5]) and ([9]) allow the Bosc- Hubbard chain to be cre- 
ated in a single resonator. It is interesting to consider 
what other many-body systems might be engineered in 
a similar way, and whether one could use an optical 
rather than a mechanical resonator to hold the result- 
ing lattice. To address the latter question first, we note 
that the asymmetry between the mechanical and opti- 
cal/superconducting resonators is that given by the in- 
teraction in Eq.([2]). All we have to do to generate a 
Kerr non-linearity for the cavity is to adiabatically elimi- 
nate the nano-resonator. This is the basis of the schemes 
in [29|, [30[. However, to form a lattice we need to selec- 
tively couple only one nano-mode to each optical mode, 
and do this for a number of optical modes. This would 
be relatively simple if one could modulate the coupling 
strength g n , which might be possible,for example, using 
a membrane-in-the-middle scenario [31| . Once we have 
created a Kerr Hamiltonian for each mode, we can se- 
lectively couple the modes together using the technique 
described above for nano-resonators. 

The main limitation in constructing lattice Hamiltoni- 
ans stems from the need to create higher-order nonlin- 
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earities perturbatively. Because of this, the higher the 
non-linearity the weaker its strength. Nevertheless, the 
third-order (x' 3 -*) nonlinearities we have considered here 
do give further interesting possibilities. As shown in poj ]. 
the cross-Kerr coupling, a^abfy, between two modes a 
and b can be engineered in precisely the same way as 
the Kerr for a single mode, by coupling two nano-modes 
to the same pair of optical modes. Another intriguing 
possibility is to use the three-way coupling between two 
super-modes and the nano-mode in Eq. to simulate a 
lattice with three-body interactions. Finally, one can use 
the on-site Kerr and cross-Kerr nonlinearities to simulate 
spin lattices in the Boltzmann (weak-coupling) regime. 
Here the cross-Kerr gives the a z <g> a z spin-spin interac- 
tion, and the linear coupling gives the appropriate weak- 
coupling (RWA) version of a x ® a x . This would allow one 
to create both integrable and non-integrable spin lattices, 
and investigate the behavior of thermalization and typ- 
icality as a function of lattice size [Hf, as well as the 



intcgrability/non-integrability transition [33 1. 

The configuration we have described above can also be 
used for quantum computing where the modes of a single 
resonator are the qubits or qdits. An advantage of this 
scenario is the ability to couple any qubit to any other 
qubit via the driving fields. The Kerr non-linearity, if 
sufficiently strong, allows the lowest two Fock-states to 
be used as a single qubit, by separating them from the 
higher states. In this case all singlc-qubit gates can be 
obtained by linear driving and the Kerr nonlincarity, and 
the cross-Kerr provides a controlled phase gate, complet- 
ing the universal set [34[. Alternatively, the Kerr non- 
lincarity, combined with the linear mode coupling and the 
ability to squeeze (achieved by parametric driving (35[) 
is universal for continuous variable computing [36]. 
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